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ABSTRACT 

We present an integral field spectroscopic study of two nearby Luminous Infrared Galaxies (LIRGs) that 
exhibit evidence of widespread shock excitation induced by ongoing merger activity, IC 1623 and NGC 3256. 
We show the importance of carefully separating excitation due to shocks vs. excitation by HII regions and the 
usefulness of IFU data in interpreting the complex processes in LIRGs. Our analysis focuses primarily on the 
emission line gas which is extensive in both systems and is a result of the abundant ongoing star formation as 
well as widespread LINER-like excitation from shocks. We use emission-line ratio maps, line kinematics, line- 
ratio diagnostics and new models as methods for distinguishing and analyzing shocked gas in these systems. 
We discuss how our results inform the merger sequence associated with local U/LIRGs and the impact that 
widespread shock excitation has on the interpretation of emission-line spectra and derived quantities of both 
local and high-redshift galaxies. 
Subject headings: Astronomy 



1. INTRODUCTION 

Understanding the physical processes involved in the for- 
mation and evolution of Galaxies remains a key problem in 
modern astronophysics. As observers strive to interpret data 
on the first generations of galaxies and theoreticians attempt 
to recreate the assembly and evolution of galaxies, studies of 
the nearby universe are needed to constrain key physical pro- 
cesses and global parameters that inform both observations 
and models of galaxies in the distant Universe. In this regard, 
Luminous Infrared Galaxies (LIRGs) and Ultraluminous In- 
frared Galaxies (ULIRGs) provide important observational 
constraints to models of distant galaxies. Although rare lo- 
cally, IR-luminous galaxies and the merger activity associated 
with this class of systems become increasingly important at 
high r edshift, contributing (ElbaFetaT|[2002; Le Floc'h et"aT] 

2005) . Not only does the improved spatial resolution allow us 
to study the relative importance of physical processes in the 
various components of the U/LIRG system, but also we can 
study the physics over the full array of merger states. Many 
recent studies across multiple wavelengths are beginning to 
focus on resolving the detailed ionization structure, kinemat- 
ics and power sources of local U/LIRGs in order to further 
our understanding of galaxy evolution both locally and at high 
redshift. 

The local population of LIRGs (which range from 10 11 
L Q <L/a < 10 12 L ) have been found to represent a va- 
riety of objects from particularly dusty isolated star-forming 
galaxies to progressed mergers or late-tage mergers. Most fre- 
quently, however, U/LIRGs are often associated with merger 
activity which drives the intense circum nuclear star formation 
that fuels the powerful IR luminosity (IVeilleux et al.l 11995b 
Scoville et aljioOOHArribas et al.l2004tlAlonso-Herrero et alJ 

2006) . The gas in the outskirts of each system is torqued 
into the center which drives the ongoing star for mation as 
the merger continues, dBarnes & H ernquist 1996). As the 



galaxies move along the Toomre & Toomre (1972) merger se- 
quence the IR luminosity is generally found to i ncrease along 
with galactic wind activity (Rupk e et alJ 12005b . At the ter- 
minal phase of the merger sequence galaxies have in general 
become ULIRGs, and appear as a coalesced or nearly coa- 
lesced system. Finally, these become dominated by a pow- 
erful AGN as their star formation becom es quenched as the 
available gas res ervoir becomes depleted dGenzel et al.lll 99 8k 
lYuan et alJl2010h . 

As this merger sequence progresses, tidally induced gas 
motions and outflows from galactic winds become increas- 
ingly common. These forces can have an influence on the 
emission line gas throug h shocks that are induced by the 
large- sc ale gas flows (lArmus et all 11989b iHeckman et al.l 

1996b iColina et all 12005b 
excitation may contami- 
nate line ratios used to determine metallicity, star formation 
rate and power source in a galaxy. Depending on the 
geometry and kinematics of the gas and dust the ionization 
of the observed emission line gas can be d ominated by 
shock excitation induced bv either tidal flows (|Cplina et al.1 
12005b iFarage et al.1 1201 Ob iMonreal-Ibero et al.1 120101) or 
galac ti c winds dVeilleux & Rupkd |2002b IVeilleux et alJ 
2003b iLiparietal l l2004at ISharp & Bland-Hawthornl 12010b 
Rich et all 12010b iMoiseev et alJ 120101) . Early work using 
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narrow band or tunable filters and Fabry-Perot instruments 
proved useful for invest igating these processes ove r large 
swaths of a galaxy (IShopbell & Bland-Hawthornl 11998b 
IVeilleux & Rupkell2002t lCalzetti et alj|2004l) . However, a de- 
tailed characterization of the shock excitation and its effects 
on emission line spectra is more complex, necessitating the 
use of integral field spectroscopy (IFS). 

The advent of integral field units (IFUs) has led to a rapid 
multiplication of observational evi dence for extended shock 
excitat ion in nearby galaxies. [Sharp & Bland-Hawthornl 
(2010) conducted a study of several galactic winds and found 
widespread shock-excitati on in the extended emission asso- 
ciated with the outflow. IMonreal-Ibero et akl ((2006) iden- 
tified extended LINER-like emission resulting from shocks 
in a small IFU sample of ULIRGs and followed up with a 
more thorough IFU study of tidally induced shock excitation 
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FIG. 1. — Archival HST ACS data of IC 1623 and NGC 3256 with our WiFeS pointings overlaid and 2 kpc bar for scale. The images shown are F435W + 
F814W, effectively Johnson B+I images. This filter combination distinguishes the intense star formation from the prominent dust lanes which in both galaxies 
heavily obscure the secondary nucleus of the merging pair. 



in nea rby LIRGs dMonreal-Ibero et al. 20101) . iFarage et ail 
(2010) discovered extended shocks caused by ga s accreting 
onto a giant brightest cluster galaxy (BCG). Finallv lRich et alj 
(2010) found extended shock excitation caused b y a galactic 
wind in the M82-like galaxy NGC 839. In both lFarage et al l 
(2010) and lRich et afl OQlO). new slow shock models were 
employed to analyze the shocked gas. In all of the above cases 
shock excitation exhibits characteristics of extended LINER - 
like emission. 

In this paper we present a detailed analysis of the 
widespread shock excitation in two galaxies drawn from a 
larger IFU sample of nearby U/LIRGs: the Wide Field Spec- 
trograph (WiFeS) Great Observatory All-Sky LIRG Survey 
(GOALS) sample. We discuss these two systems, IC 1623 
and NGC 3256 in § [2] and our observations, data reduction 
and line-fitting in § [3] We present the properties of the emis- 
sion line gas fo r these two systems in § |4] including line ratio 
maps in § 14. II and emission line diagnostic diagrams in 34.21 
We discuss the distribution of velocity dispersions for each 
system in § [5] We combine the observed quantities and use 
them to separate the shocked gas from the HII region gas and 
provide an analysis of our results in § |6] In § [7] we discuss 
the implications of our results for the interpretation of emis- 
sion line spectra and investigate the effect that extended shock 
excitation can have on quantities derived from single aperture 
spectra of nearby and high-redshift galaxies. Finally we give 
our conclusions in § [8] 

Throughout this paper we adopt the cosmological parame- 
ters H =70 km s^Mpc" 1 , f2 y=0.72, and £^=0.28 , based on 
the five-year WMAP resu lts Hins haw et al.l (120091) and con- 
sistent with the lArmus et al.1 (12009) summary of the GOALS 
sample. Systemic heliocentric v elocities adjusted a djusted us- 
ing th e 3-attractor flow m odel of M ould et al.l (120001) are taken 
from lArmus et al.l (120091) . 

2. SAMPLE 

The systems analyzed in this paper are drawn from a larger 
integral field spectroscopic survey of the GOALS sample. 
GOALS is a multi-wavelength survey of the brightest 60/im 



TABLE 1 



Galaxy 


z 


RA. 


Decl. 


B Mag. 


log(L/ fl /L Q ) 


IC 1623 


0.0201 


01:07:47 


-17:30:27 


15.04 


11.71 


NGC 3256 


0.00935 


10:27:51 


-43:54:14 


11.83 


11.64 



extragalactic sources in the local universe (log(Lm/ X©) > 
11.0), with redshifts z < 0.088 and is a complet e subset of the 
IRAS Revised Bright Galaxy Sample (RBGS) (ISanders et alJ 
2003). Objects in GOALS' cover the full range of nuclear 
spectral types and interaction stages and provide excellent 
nearby laboratories for the study of galaxy evolution. 

We chose IC 1623 and NGC 3256 for a more detailed anal- 
ysis since they both show evidence of widespread shock ex- 
citation. Table 1 shows the global properties of the two sys- 
tems, both of which are ongoing mergers of two massive spi- 
ral galaxies. Both systems are fairly luminous LIRGs and 
neither shows evidence of an embedded AGN. Archival HST 
data of the systems is shown in fig. [TJ Both systems show a 
complex optical morphology and contain one galaxy nucleus 
that is heavily extinguished by dust lanes. Both galaxies are 
classified as close mergers (projected separation between the 
merging nuclei of l ess tha n 10 kpc) when using the scheme of 
Veilleux & Rupke (2002), which utilizes classifications based 
on numeric al simulations of mergers (Ba rnes & Her nquist 
[T99l[T996h . 

IC 1623: IC 1623 (also commonly referred to as VV 114 
and with the IRAS designation IRASF 01053-1746) consists 
of two galaxies: IC 1623 A in the west and the more heavily 
obscured IC 1623 B in the east. The dust lane seen in figureQ] 
covers the eastern nucleus, w hich is actually brighter in the 
near-IR than the western d isk (Kno p et aill 19941 : IDo von et alj 
119951: iScoville et al.ll2000l) . The galaxy nuclei are separated 
by approxim ately 5 kpc, firmly classify ing the pair as a close 
merger in the Veilleux & Rupke (2002) classification scheme. 
Interestingly, the heavily obscured eastern galaxy has two 
very bright sources of similar magnitude in the near-IR, rais- 
ing the possibility that IC 1623 B itself is a merger. However, 
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it may be that there is simply a second region of heavily ob- 
scured star formation unassociated with the true nucleus. Be- 
cause of the extensive dust lane, the majority of the optical 
emission in our observations is dominated by the star forma- 
tion in the western galaxy. Neither of the gala xies exhibits 
evidence of an AGN, ob scured or otherwise (Dov on et al.l 
119951; IVeilleux et all 119951) . Extended tidal featur es can be 
seen in optical and near-IR imaging, though the low surface 
brightness of these features makes observations of them pro- 
hibitively time-consuming. The galaxy pair is at a redshift of 
z— 0.0201, corresponding to a 1 "projected distance of —385 
pc. 

NGC 3256: NGC 3256 (IRAS designation IRASF10257- 
4339) is a late-stage merger with nuclei separated by ap- 
proximately 1 kpc, though it has not yet formed a single 
nucleus and the r efore cannot be classified as a post-merger. 
ISakamoto et alJ (120061) n ote that the long tidal tails seen 
in op ti cal and HI images ( de Vaucouleurs & de Vaucouleursl 
I196U lEnglish et al.l 120031) suggest "a pro grade-pro grade 
merger of two gas-rich spiral galax ies of similar size" 
dToomre & Toom re 1972; IWhitel 11979b . The bright face-on 
galaxy seen in figure [T]has a separated, heavily obscured sec- 
ondary nucleus seen in x-ray, IR and radio lying —5" south 
of the bright optical nucleu s (iBoeker et al.l 1 1 9971: iNeff et aD 
l200l ISakamoto et ail I2006T) . Neff et alJ d2003lT posited that 
the two nuclei may be low-luminosity AGN based on the ra- 
dio to x-ray ratios, though they also note that a large number 
of supernova remnants could also accoun t for the observed ra- 
tio. However subsequent observations by lLipari et al.l (12004b) 
of galactic bubbles associated with previous supernovae in 
the core of NGC 3256 support the latter case and rule out 
any significant contribution from even a low-luminosity AGN. 
This supernova activity in the core of NGC 3256 is also 
driving a galactic wind measured in neutral absorption and 
ionized emission previously associated with shock heating 
dArmus et al.|[l989t iHeckman et al.lfl990l l2()0"ot Opari et all 
2000). NGC 3256 lies closer than IC 1623 at z~ 0.0094, 
which means that our data probe a finer spatial scale of 
1"- 185 pc. 

3 . OB SERVATIONS & DATA REDUCTION 

Our data were taken using the Wide Field Spectrograph 
(WiFeS) at the Mount Stromlo and Siding Spring Observatory 
2.3 m telescope. WiFeS is a dual beam, image-slicing integral 
field spectr ograph commission ed in May 2009 and des cribed 
in detail by Do pitaetalJ (120071) and lDopita et all d2010t) . Our 
data consists of separate blue and red spectra with wavelength 
coverage of -3500-5800 A & -5500-7000 A respectively 
and a resolution of R = 3000 (blue) and R = 7000 (red), cor- 
responding to a velocity resolution of 100km s _1 at H/3and 
40km s _1 at Ha . Observations of IC 1623 were carried out 
on August 20, 2009. Two separate pointings were observed 
for 40 minutes in total for each pointing. NGC 3256 was ob- 
served on March 15, 2010 with two individual pointings for 
60 minutes each. Each galaxy was observed using two point- 
ings as shown in Q] with a large enough overlap region to 
facilitate mosaicing. The effective field of view of the data in 
both cases is thus approximately 45" x 38", providing well 
over 1,000 spectra in each set of data cubes. 

The data were reduced an d flux calibrat e d usin g the WiFeS 
pipeline, briefly described in Dopi ta etldl (12010). which uses 
IRAF routines adapted primarily from the Gemini NIFS data 
reduction package. A single WiFeS observation consists of 25 
slit spectra (and an additional 25 contemporaneous sky spec- 



tra if the data was taken in nod and shuffle mode). Each in- 
dividual observation was reduced into a data cube using the 
process described below. 

Bias subtraction is slightly complicated due to temporal 
changes in the quad-readout bias levels and a slope in the bias 
level on each WiFeS chip. Bias frames are taken immediately 
before or after each observation or set of observations and a 
2-D fit of the surface of this bias frame is subtracted from the 
temporally nearest object data in order to avoid adding ad- 
ditional noise to the data. Any resulting residual in the bias 
level is accounted for with a fit to unexposed regions of the 
detector. 

Quartz lamp flats are used to remove pixel to pixel vari- 
ations from the response curve of the chip and twilight sky 
flats are used to correct for illumination variation along each 
of the slitlets. Spatial calibration is carried out by placing a 
thin wire in the filter wheel and illuminating the slitlet array 
with a continuum lamp. This procedure defines the center of 
each slitlet. The individual spectra have no spatial distortion 
because the camera corrects the small amount of distortion in- 
troduced by the spectrograph. Thus only low-order mapping 
of the slitlets is required. 

CuAr and NeAr arc lamps are used to wavelength calibrate 
the blue and red spectra respectively. Arc lamp data were 
taken throughout each night of observing to account for any 
change in the wavelength solution. Each of the 25 slitlets is 
then rectified by the pipeline into a full data cube (one for each 
arm) and sampled to a common wavelength scale for each tar- 
get. The spectral resolution achieved within the red data cubes 
is 46 ± 6 km s -1 . We use this value to remove instrumental 
broadening from the derived velocity dispersions. 

Telluric absorption features were removed from the result- 
ing red data cubes using observations of B-stars or featureless 
white dwarfs taken at similar air mass. The effects of differ- 
ential atmospheric refraction are calculated and corrected by 
the pipeline for the blue data cubes. 

Flux calibration of each individual data cube was carried 
out using flux standards observed throughout each night. The 
white dwarf EG 131 was used to flux calibrate IC 1 623 and the 
white dwarf L 745-46a was used to flux calibrate NGC 3256 
dBesselllll999h . 

Finally, the individual reduced, flux-calibrated data cubes 
are median combined and sampled to a common spatial grid 
using overlapping features found in each pointing. Observa- 
tions of NGC 3256 were binned on-chip by 2 pixels in the 
spatial direction in order to increase signal to noise and pro- 
duce square spatial elements 1" x 1". Observations of IC 1623 
were binned after reduction for the same reasons. This means 
that for both systems the scale of the data cubes and of the 
resulting maps presented in this paper are l'Vspaxel. 

Combined data cubes were aligned astrometrically by com- 
paring a pseudo 'r-band' image generated using the red spec- 
trum from each spaxel with HST images from the archive. 

3.1. Spectral Fitting 

We analyzed every spectrum using an automated fitting rou- 
tine written in IDL, U HSPECFIT, which is based on the code 
created by Uabran Zahid et al.l (12010 1 ) and is also employed 
by and described in iRupke et al.l ( 1201 Obi) and iRich et all 
(2010). Our routine fits and subtracts a stellar continuum 
from each spectrum using pop ulation synthesis models from 
iGonzalez Delgado et al.l d2005[) and an IDL routine which fits 
a linear combination of stellar templates to a galax y spectrum 
using the method of Moustakas & Kennicutt] ( 120061) . 
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FIG. 2. — Example of emission line fitting routine results for a single spectrum from IC 1623. The upper panels show individual emission lines of interest with 
data in black, total fit in red, continuum in green, and residuals in blue. The bottom panel shows a zoomed in region around [Nil] and Ha with each Gaussian 
component plotted as a dotted line and the total fit plotted as a solid line. Wavelengths are plotted in the rest frame with respect to the narrow component. 
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After the continuum is removed, lines in the resulting emis- 
sion spectra are fit using a one or two-component Gaussian, 
depending on the goodness of fit determined by the routine. 
Resulting fits were inspected to ensure the fitting code did not 
fail. All emission lines are fit simultaneously using the same 
Gaussian component or components. Both continuum and 
emission lines were fit using the MPFIT package, which per- 
forms a lea st-squares an alysis using the Levenberg-Marquardt 
algorithm (Markwardt 2009). Example fits are provided in 

fig-El 

The errors in the parameters used in fitting the emission 
lines are calculated by the fitting code. This includes the 
gaussian widths (velocity dispersions) and ratio of dispersions 
between the blue and red spectra, individual gaussian peaks, 
redshift and any slight deviation between the blue and red 
wavelength calibrations. Emission line fluxes used in derived 
quantities are subject to a cut of minimum signal-to-noise of 
5, which is calculated by comparing the peak of the Gaus- 
sian to the noise in the continuum. The lower resolution of 
the blue spectra contributes to the uncertainty of the gaussian 
parameters-for instance in some cases 1 or 2 emission line 
components for a blue emission line may provide a similar 
goodness of fit. By minimizing the \ 2 for a ll of the lines si- 
multaneously we increase slightly the overall uncertainty in 
the velocity dispersion while maintaining the ability to ana- 
lyze multiple components fits. 

In many cases it is crucial that more than one Gaussian is 
used in the emission line fit. Our analysis relies on a careful 
decomposition of multiple velocity profiles, which are com- 
mon in the kinematically complex family of U/LIRGs. An 
example can be seen in fig. [2] where two distinct compo- 
nents have been fit. In the case shown, there is an obvious 
red-shifted component with a broader profile than the primary 
component. The line ratios between the two components are 
noticeably different as well. 

4. EMISSION LINE GAS PROPERTIES 

Emission line ratios have proven a useful tool for probing 
the sources of ionizing radiation in galaxies. The primary 
factors influencing emission line strengths are the shape of 
the extreme ultraviole t (EUV) spectrum and the abundance o f 
the emission line gas dKewlev et alJl200ll : lDopita et alj|2000l) . 
Thus with IFU data, a line ratio maps and diagnostic dia- 
grams can be created as a useful tool in interpreting the power 
sources not just in the nucleus of a galaxy but also across the 
galaxy as a whole. 

4.1. Line Ratio Maps 

In figure 3, we examine maps of the ratios of the total flux 
of [Nil] 6583 A, [SII] 6717,6731 A, and [OI] 6300A to Ha 
which are sensitive to metallicity as well as the ionization pa- 
rameter as a first step towards understanding the processes at 
work in IC 1623 and NGC 3256. 

IC 1623: Most of the optical emission line flux in IC 1623 
is dominated by the eastern galaxy (IC 1623 A), which is less 
extinguished by dust. The emission line ratios in this por- 
tion of the merging system are lower than in the rest of the 
system, tracing the star formation evidenced by the less ex- 
tinguished HII regions and clusters seen in the archival HST 
images. There is no evidence of a metallicity gradient in the 
[NII]/Ha ratio map, which is in agreement with recent work 
on the flattening of met all icity gradients in me rging systems 
bv iRupke etalT ( l2010allb1) : lKewlev et alj feOlCft . 



Moving outward from the main regions of unobscured star 
formation in IC 1 623 A we see an overall rise in all line ratios . 
Towards the northern and southern parts of IC1623 A, the line 
ratios increase to levels inconsistent with excitation by HII re- 
gions alone until the signal to noise in emission is too low. 
Especially striking are the high line ratios correlated with the 
dust lane east of IC 1623 A which obscures the better part of 
IC 1623 B. The emission line ratios show consistently high 
values across this dust lane, especially in [OI]/ Ha , which is 
a particularly good tr acer of shock excitation (Fara ge et alJ 
l2010tlRichetal.ll2010l) . The total [OI] flux in this region ap- 
proaches 1/4 to 1/3 the total Ha flux, while the total [Nil] flux 
is as much as 1/2 to 2/3 the total Ha flux. 

NGC 3256: The line ratio maps of NGC 3256 show simi- 
lar behavior compared to IC 1623. Again the areas of unob- 
scured star formation correspond to lower line ratios, tracing 
the ongoing star formation as the dominant source of ionizing 
photons. A key difference is that NGC 3256 shows a region 
of increasing line ratios in a small central part of the region 
of low line ratios. This increase correspon ds to the previ- 
ously observed galactic wind in NGC 3256 (Heckm an et al.l 
1990) and our spectra in this region exhibit strong, broad blue- 
shifted NaD profiles (fig. |5). The increasing emission line 
ratios in this case are due to the additional presence of shock 
excitation from the galactic wind. 

Outside of both the central region influenced by the galac- 
tic wind and the primarily star formation dominated por- 
tions of the galaxy the line ratios again tend to increase, 
with [NII]/Ha reaching nearly unity and correspondingly high 
[SII]/Ha and [OI]/Ha in the outer spatial portions of our data 
cubes. These regions of strong excitation correlate strongly 
with dust lanes evident in the HST archival data and are quite 
distant from the core of the galaxy. 

4.2. Diagnostic diagrams 

In order to better understand the power sources at play, 
we turn to line ratio diagnostic diagrams. Diagnos- 
tics using [Nil]/ Ha,[SII]/Haor [OI] /Ha ratios against the 
rOI II]/H/3 ratio were first employe d by iBaldwin et ail (11981b 
and IVeilleux & Osterbrockl (11987b to distinguish the source 
of ionizing photons in global or nuclear galaxy spectra. This 
diagnostic scheme has been successfully applied to Integral 
field spectra (as well as narrow-band and tunable filter data) 
to investigate emission line gas over entire galaxies or regions 
of a galaxy to probe the distribution of ionizing sources (e.g. 
Shopbell & Bland-Hawth ornll998HSharp & Bland-Hawthornl 
2010tlRich et alj|2010t) . To separate the ionizing sources with 
IFU data, individual spectra from each spaxel are plotted on 
the traditional diagnostic diagrams. In this way the ensemble 
of points provides information about the various processes at 
work in different portions of the galaxy rather than a single 
global measure. 

We plot diagnostic diagrams in figure |4] and include all 
emission line components from each spaxel with S/N> 5 in 
all 6 diagnostic lines (H/3 , Ha , [OIII] 5007 A, [OI] 6300A, 
[Nil] 6583A, [SII] 6717,673lA). In all diagrams, the lower 
left-hand portion of the plot traces photionization by HII re- 
gions. The solid curved line traces the upper th eoretical limit 
to pur e HII region contribution measured by iKewley et al.l 
d2001l) and the dotted line in the [Nil] diagnostic provides 
an empirical upper limit to the HII regio n seque nce of SDSS 
galaxies measured by Kauff mann et alj d2003l) . The region 
lying between these two lines represents objects with a com- 
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FIG. 3. — From left to right, emission line ratio maps of [NII]/Ho , [SII]/Ha and [OI]/Ho for IC 1623 (top) and NGC 3256 (bottom). Ha contours have been 
overlaid, and the peak in Ha flux has been marked with a small red "x". These maps are a sum of the total emission line flux in both profiles and trace the 
dominant ionization mechanism in different portions of the galaxy. Generally speaking, the stronger the line ratio the further the departure from an ionizing 
spectrum purely due to star formation. The regions most dominated by ongoing star formation appear bluer in the above maps. Maps are oriented with respect to 
the original observations (cf. fig. 9) 



Line Ratios for single and double components 



02. 
X 



o 



02. 
X 



O 



1.5 

1.0 

0.5 

0.0 

-0.5 

-1.0 
1.5 

1.0 

0.5 

0.0 

-0.5 

-1.0 



IC1623 

' ><N AGN 

hii * .•••.*A 

■ ■ ■ ■ 1 1 1 ■ ■ i ■ ■ ■ ■ i ■ ■ i ■ i jIIS .ii. 


— ■ — i — ■ — i — i — ■ — i — i — i — ■ — ■ — i — i — i — ■ — ■ — i — ■ — 

Sy2 

"- . » / 

Hl1 \ LINER . 
■ 1 ■ ■ i ■ 1 i ■ ■ ■ i' ■ i ■ ■ 1 ■ 


Sy2 

<» v •• 

Hl1 \ LINER 
■ 1 ■ ■ ■ ■ 1 ■ i ■ ■ 1 ' i i ■ 1 ■ ■ i ■ 1 ■ 


NGC 3256 




X ✓ 
•*\ ij* 

• ' ' \ ' 

• 1 



-2.0-1.5-1.0-0.5 0.0 0.5 
log([NII]/Ha) 



1.0 -0.5 0.0 0.5 -2.0 -1.5 -1.0 -0.5 0.0 
log([SII]/Ha) log([OI]/Ha) 



FIG. 4. — Emission line ratio diagnostic diagrams of each individual fitted emission line component for each system. Boundaries are taken from Kew lev et al.l 
1 2006). We apply a minimum S/N cut of 5 in every diagnostic line. The resulting points in each diagram reflect a mix of ionization by HII region emission and 
by shock excitation. The points lying in the HII region portion of the diagram also reflect the abundance of the gas: the HII region gas in NGC 3256 is higher 
metallicity resulting in a different overall shape to the distribution of points in the diagnostic diagram. 



7 



Na D + He I 

2.0x 1 0" " I 1 1 1 1 1 1 r 




Ol I I I I I I I 

5928 594J 5958 ^ J^j 5928 5943 5958 



Ha + [Nil] 















TV J 


A : 



6609 6626 6642 X Q ^ S (A) 6609 6626 6642 



FIG. 5. — Line profiles from 3" extracted apertures from NGC 3256 for 
the brightest optical nucleus (left) and an off-nuclear region dominated by 
star formation (right). The nuclear spectrum shows a large blue-shifted Na D 
component as well as a corresponding blue tail in emission both caused by 
the large galactic wind driven by the nuclear starburst. 

posite spectrum mixing HII region emission with a stronger 
ionizing source. Contribution from a Seyfert AGN pushes 
points into the upper region of the diagrams while more 
LINER-like emission lies to the right hand side of the dia- 
grams, as shown by the dividi ng line in the [Sill a nd [OI] 
diagnostics, as established in Kewley et alj (l2006h . New 
lower velocity shock models of order 100-200 km s _1 trace 
line ratios from the upper edge of HH-lik e emission into the 
LINER region of the diagnostic diagrams (Farage_etalJ[2010; 
iRich et al.ll20l0h . 

5. VELOCITY DISPERSIONS 

Where more than one component is fit to emission lines 
ordering the components over a whole IFU data cube is non- 
trivial. The difference in velocity dispersion and flux where 
two components are fit varies depending on the conditions in 
the region of the galaxy observed and is dependent on several 
factors including the relative flux, source of ionization and 
extinction within each component. In some cases two compo- 
nents of comparable flux and velocity dispersion are observed, 
while in other regions two components of significantly differ- 
ent flux, flux ratios and velocity dispersion are seen (cf |2). In 
this section we focus on the velocity dispersion of individual 
emission line components as the most useful way of probing 
the multiple components observed 

If shock excitation is important in these regions showing en- 
hanced [OI]/Ha, [NII]/Haand [SII]/Haratios, then we would 
expect to see the presence of these shocks reflected in in- 
creased velocity dispersions in these areas. As the emission 
line ratios increase in a manner consistent with shocked gas, 
the velocity dispersion should increase to velocities consis- 
tent with the speed of the shocks, which are of order 100-200 
km s _1 . Lower emission line ratios consistent with gas in HII 
regions should exhibit lower velocity disp ersions, typically of 
a few 10s of km s ~ 1 ( Ep inat et al.l 120 1 Oh . In a similar study 



of LIRGs with IFUs, iMonreal-Ibero et al.l d2006l l2010b ob- 
served the correlation between excitation and velocity disper- 
sion consistent with widespread shock excitation in a larger 
sample of nearby U/LIRGs and associated this excitation with 
tidally induced shocks. 
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FIG. 6. — Distribution of emission line velocity dispersions from every in- 
dividual Gaussian component fit. The bimodal distribution of velocity dis- 
persions between the low velocity (unresolved HII region) and high velocity 
(shock-excited) components can be seen quite clearly and the cutoffs between 
low and high velocity dispersion are marked with dotted lines. The bimodal 
distribution in NGC 3256 is less sharply defined due to the combination of 
outflow and inflow. 

We are able to more accurately separate the lower veloc- 
ity dispersion (low-c) emission consistent with HII regions 
and the higher velocity dispersion (high-cr) with our multiple 
component emission line fitting, as large portions of the gas 
in both systems exhibit a composite of activity both in emis- 
sion line ratios and in velocity dispersion. For instance, in IC 
1623 the gas falling on to the dust lane is well fit by a single 
fairly broad component, while gas elsewhere in the galaxy is 
mostly fit by either two narrow components or a narrow and 
broad component consistent with the complex velocity field 
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Line Rotio vs. Velocity Dispersion 
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FIG. 7. — Plot of line ratio vs. velocity dispersion for each individual emission line component. Regions of emission line gas dominated by shock excitation 
show both higher emission line ratios and higher velocity dispersion resulting in the correlation seen above. A sequence of narrow component, low-excitation 
points corresponds to HII region gas of varying ionization parameters and metallicities. There is less scatter in the [NII]/Hq line ratio owing to the higher signal 
to noise of that line, though the broader spread of [SH]/Hct and especially [OI]/Ha makes these ratios a better tracer of shock excitation. 



and multiple-phase gas in this merger. 

In fig.[6]we plot a histogram of velocity dispersions for each 
system for every component fit by our routine. The spectrum 
from every spaxel was fit with either one or two Gaussian 
profiles with unique widths, and all of those components are 
plotted here. We only plot values from spaxels with emission 
lines of S/N> 5 in all of our diagnostic lines, H/3, [OIII], 
[OI], [Nil], [SII] and Ha. We have removed the effect of 
instrumental broadening from our velocity dispersions. 

In both systems there is a clear peak at low-cr (~ 40 km s _1 
consistent with a large contribution from HII region emission 
due to the substantial amount of ongoing star formation. In 
the case of IC 1623 there is a clear bimodality between the 
low-cr (HII region) and high-cr (shocked) components. The bi- 
modality in NGC 3256 is still present but appears less sharply 
delineated owing to a combination of shock excitation by the 
galactic wind and by bulk gas motions due to the ongoing 
merger. 

We use the histograms in fig.|6]to establish a cutoff between 
the low-cr and the high-cr gas. For each system dashed lines 
show the cuts in fig. [6] which separate the emission line com- 
ponents into low-cr, transition and high-cr to better trace the 
different power sources. For IC 1623, we consider all gas with 
cr < 60 km s _1 as low-cr and all gas with a > 90 km s _1 as 
high-cr with the transition falling in between these two values. 
Similarly, we establish cutoffs of low-cr at < 50 km s _1 and 
high-cr at > 65km s^forNGC 3256. 

Figure [7] shows velocity dispersion vs. emission line ra- 
tio for each individual component. Correlation between 
these two measures of shock excitation was explored us- 
ing single-compone n t fits t o IFU data of several LIRGs by 
Monreal-I bero et aTl (120101) . The increasing velocity dis- 



persion correlated with increasing emission line ratio corre- 
sponds to the dominance of shock excitation at high-cr, while 
the narrow sequence of low-cr points (bottoming out at the 
instrumental resolution) corresponds to emission from HII re- 
gion gas. 

6. ANALYSIS 

In figure |8] we color-code the points in our diagnostic dia- 
grams from figure |4] using the low-cr and high-cr cutoffs for 
each system in order to separate the HII region gas from 
the shocked gas. The correlation of velocity dispersion with 
line ratio strength is still clearly seen, most obviously in 
the [Nil] diagnostic owing to the higher S/N of that line. 
More importantly, there is a clear sequence from Hll-region 
to LINER in the [SII]and [OI]diagrams with a correspond- 
ing change from low-cr to high-cr gas tracing the influence 
of shocks. The locus of the HII region emission-line gas 
points in the [Nil] diagnostic diagram begins in the mid- 
dle/low metallicity region for IC 1623 and in the higher metal- 
licity region for NGC 3256, corresponding to the abundances 
of the HII regions in those galaxies. 

Color coding the points accentuates the mixing sequence 
from pure HII region emission to pure shock excitation al- 
ready seen in the emission line diagnostic diagrams. By com- 
bining emission line diagnostic ratios and velocity dispersion 
analysis of individual emission line components we can per- 
form a fairly conservative separation of the HII region gas and 
the shocked gas of each system and analyze them separately. 

In fig. [9] we plot the line of sight velocity of the low-cr and 
the high-cr separately for both systems. The first thing of note 
is the presence of a high-cr component over nearly the entirety 
of both systems-clearly shocks are widespread in these galax- 
ies. We have excluded some points from the central regions 
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FIG. 8. — Same plot as fig. [4] but with points color coded based on velocity dispersion cutoffs established in section [5] Low-cr points are plotted in blue, 
transition points in yellow and high-cr points in red. Clearly the low velocity (HII region) and the high velocity (shocked) components are characterized by 
different ranges of line ratios. 

km s _1 and a broader narrow component such that both com- 
ponents have strong line ratios and are considered high-cr and 
thus excluded from the velocity plots. Finally, there is a region 
~ 8"in diameter of high-cr component strongly blue shifted- 
up to nearly 200 km s -1 -lying northeast of the nucleus in a 
region with little dust and several closely-packed star clusters, 
possibly a second localized outflow. 



of each system where our fitting routine was unable to satis- 
factorily fit 2 Gaussians to the complex line profiles or where 
both emission line components are broad enough to exceed 
the shock cutoff. In both systems the low-cr gas traces the 
HII regions which follow the rotational motion of the parent 
galaxies. 

In IC 1623 the low-cr gas is primarily dominated by the 
rotation of the eastern galaxy and is comparable to the over- 
all shape of the rotation observed in the CO gas observed by 
lYunet al.l(fl994l) . There is a gap in low-cr gas consistent with 
the dust lane seen in the HST images, as the shock excitation 
is correlated with the dust lane. In this region our emission- 
line analysis only traces the tidally-induced shock excited gas 
between the two systems because the underlying star forma- 
tion of the western galaxy is completely extinguished by dust. 
The kinematics of the shocked gas are quite distinct, exhibit- 
ing very little correlation with the rotation seen in the low-cr 
gas. The shocked gas in the northern portion of IC 1623a is 
blue shifted up to 200 

The kinematic distinction between the HII region gas and 
the shocked gas is more subtle in NGC 3256. Where both 
a shocked component and an HII region component are both 
measured, however, the shocked gas is mostly blue-shifted. 
Over the larger part of the visibly star forming regions of 
NGC 3256 there is a mild shocked component blue shifted 
a few 10s of km s _1 . The dust lane extinguishes the HII re- 
gion gas as in IC 1623 and in spaxels where both a low-cr 
and a high-cr component are measured along the dust lane the 
shocked gas is generally redshifted a few 10s of km s _1 . In 
addition, NGC 3256 hosts a galactic wind seen q uite clearly 
in blue shifted Na D dHeckman et alJll990l l2000h . The spec- 
tra affected are within a few arcseconds of the nucleus and 
have a very broad blueshifted tail extending several hundred 



7. DISCUSSION 

Our analysis clearly separates the HII region gas from the 
shocked gas where both are detected. In both systems along 
the dust lanes and in the periphery of the galaxies where star 
formation either decreases or is extinguished the emission line 
spectra are dominated by a single, shocked component. How- 
ever, our results indicate that even in the regions where emis- 
sion is dominated by star formation there is still an underlying 
broad component associated with shocks, induced either by 
outflows generated by the circumnuclear star formation or by 
the tidal motions of the gas. 

7.1. Shock and HII region Models 

To characterize this shock excitation we employ new slow 
shock models which reproduce our me asured emission-line 
ratios. These models are introduced in Farage et al.l (12010b 
and iRich et all (120101) and will be discussed in detail in an 
upcoming paper (Kewley et al. in prep). In short we gen- 
erate slow shock models with velocities consistent with our 
observed line widths using an updated version of the Map- 
pings HI code, originally introduced in Sutherland & Dopita 
(1993). In addition, w e compare HII region models gener- 
ated using Starburst99 (Leithereret al. 1999) and Mappings 
III. Both sets of models have consistent abundances of solar 
and twice solar for IC 1623 and NGC 3256 respectively, using 
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FIG. 9. — The above figure shows the line-of-sight velocities in km s _1 of the HII region (low-cr) gas and the shocked (high-cr) gas for IC 1623 (top) and NGC 
3256 (bottom). The leftmost panel is the HST ACS image from fig.[TJfor reference. The HII region gas in both systems traces the rotational motion of the galaxy. 
The shocked gas in IC 1623 shows a dramatically different kinematic profile associated with the bulk motions of the shocked gas with components both blue 
and red shifted from the HII region component. Overall the broad component in NGC 3256 is slightly blue shifted with respect to the HII region gas in the star 
forming regions and slightly red shifted along the dust lanes. 



the abundance set of Grevess e et al.l (120101) . HII models are 
plotted with varying ionization parameter Q(H), the flux of 
ionizing photons divided by the hydrogen atomic density (the 
dimensionless ionization parameter U = Q /c), while shock 
models are shown at varying shock velocities. 

These models are plotted in figure [10] with the measured 
emission line ratios underplotted as small points. As in fig- 
ure [8] the points lying furthest to the left are more consis- 
tent with pure star formation. Our HII region models re- 
produce these line ratios. The data for IC1623 is consis- 
tent with a range in ionization parameters between 7.0 < 
log[Q{H)] < 8.0, while the higher-abundance NGC 3256 is 
characterized by a lower mean ionization parameter 6.75 < 
log[Q(H)] < 7.25 consistent with the dynamical models of 
HII regions i n different metallicity environments presented by 
Dopit a etaf] (l2006allbl) . As shocks become a more dominant 
ionizing mechanism and/or as the light from HII regions is 
extinguished by dust allowing shocks to dominate the line ra- 
tios migrate along the theoretical mixing sequence shown in 

fig. Ho] 

The shape of the mixing sequence evident in figs. [8l&[T0l 
within the diagnostic diagrams is caused by the different un- 
derlying abundances in the two systems. In IC 1623 the lower- 
abundance gas means as shock excitation begins to dominate 
the line ratios move horizontally across the diagrams to the 
LINER region. The relatively higher abundance of NGC 3256 



however creates a tilted sequence with the pure HII region 
points beginning in a different region of the diagnostic dia- 
grams. This difference is most noticeable in the [Nil] di- 
agnostic diagram where the different metallicities are most 
widely separated in line ratio space. 

Varying the fractional contribution of shock excitation 
within the parameters of our models recovers the observed 
emission line ratios quite well. It is important to note in [TO] 
that the line ratios are sensitive to as little as 20% shock- 
excitation. The observed ratios in IC 1623 covernearly the en- 
tire mixing sequence, implying that the various regions sam- 
pled in the galaxy cover a range of shock velocities, though 
most of the points near the end of the mixing sequence cor- 
respond to 100 to 140 km s _1 . In NGC 3256 the line ra- 
tios both within the mixing sequence and near the pure shock 
models are consistent with shock velocities of ~ 160 — 200 
kins" 1 . In both systems the terminus of our fastest shock 
models shows remarkable agreement with the strongest mea- 
sured emission line ratios. 

7.2. Global Energetic Considerations 

It is worthwhile to make a few simple calculations to quan- 
tify the impact that shocks have on dissipating the energy of 
the infalling gas in merging systems. Here we focus on the 
well-studied NGC 3256, where previous wide-scale observa- 
tions facilitate these calculations. 
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FIG. 10. — Here we replot measured emission line ratios (small points) with our HII region and shock models overlaid. The HII region models are larger blue 
circles connected by a black line and are plotted with varying ionization parameter (Q) and adopting abundances of solar and twice solar respectively for IC 1 623 
and NGC 3256. The fully preionized slow shock models are plotted as larger red circles and are also connected by a black line along a sequence of increasing 
shock velocity. The colored points in between show a mixing sequence from pure HII region to pure shock excitation. 



Once the shocked emission in our data has been isolated 
we can use the measured Ha flux provided solely by shocks 
to extrapolate the total shock luminosity. For this calculation 
we make the simple assumption that the total Fn a measured 
in any line components above the shock velocity dispersion 
th reshhold is cause d by shocks. We use the plot provided 
in lRich et alj (|2010) to approximate the total shock luminos- 
ity from Ln a originating in shocks. This provides an upper 
limit to the total shock luminosity of 80x Ln a ■ For NGC 3256 
this corresponds to an approximate global shock luminosity of 
Lshock ~8x 10 i2 ergs s _1 

Working with this luminosity as a starting point we next 
consider the total gas available as wel l as the mass and size 
of the system. lEnglish & Freemanl (2003); English et al. 
(2003) report an I-band (stellar) mass within 70" of the sys- 
tem of (2.5 ± 0.5) x lO lo ft, _1 M and a total HI gas mass 
of 3.5 x 10 9 /i~ 2 M Q , 75% of which resides in tails extend- 
ing to a total distance of 42ft. _1 kpc. This does not take 
into account the molecular gas mass, which may be larger. 
Adopting our value of h=0.7 provides a dynamical timescale 
of ~ 800Myr and a total binding energy for the HI gas of 
~8xl0 56 er 5 s. 

Th e star formation rate can be approximated using L h „ and 
L f p dKennicuttll 1 9981: iDopita et alj200l 120031: iCalzetti et alj 
120071 IKennicutt et alj|2009l) . iLira et alj ( 120021) estimate the 
total Fjf„ ~ 2 x 10~ 11 ergs s _1 cm -2 , with data from 
Lfp ari et al.1 ([2000) using a larger aperture than is available to 



our IFS data. Applying a total reddening correction of E(B — 
V) = 0.5, approximated using maps of the balmer decre- 
ment from our data cubes, yields Ln a ~ 1.2 x 10 43 ergs s^ 1 . 
Assuming a fractional contribution to Ha from shocks of ~ 
30%, consistent with our IFS aperture, the total star forma- 
tion rates are ~ 66M jyr and ~ 75M® jyr using Ln a and 
L/ij respectively. The addition of both extinction correction 
and an adjustment due to shock contribution yields a better 
agreement in the star formation rates. If a SFR of 75 Mq / yr, 
persists the available HI gas can fuel the starburst at this rate 
for~100Myr. 

For shocks to aid in removing kinetic energy from the in- 
falling gas, we consider the amount of time it would take for 
the shock luminosity to dissipate about twice the binding en- 
ergy of the available HI gas, or ~ lQ 56 ergs. If the entirety of 
the the shock luminosity of L s h oc k ~ 8xl0 42 ergs s _1 acts in 
this way, the kinetic energy dissipation time scale is ~ 6 Myr, 
significantly less than the dynamical time scale, implying that 
shocks may act as a very efficient mechanism for removing 
kinetic energy fr om the infalling gas , allow ing it to feed the 
central starburst. English & Freemanl (120031) . however, do not 
consider inclination effects in their mass calculations, which 
would effectively increase this timescale. 

7.3. Star Formation Rate 

The Ha flux is often used as a direct proxy for star for- 
mation ra t e in g alaxies, typically assuming the relation of 
IKennicutt! ( 119981) . If one assumes that all of the high-tr 
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shocked emission in each system is coming from pure shock 
excitation, a simple total fractional contribution from shocks 
and thus the effect on the measure of star formation rate can 
be estimated. 

Internal extinction will also impact the measurement of 
Ha -derived SFR and may affect regions of shock excita- 
tion differently than regions dominated by HII photoioniza- 
tion. When comparing the balmer decrements from our data 
for both IC 1623 and NGC 3256, the average difference be- 
tween these regions is negligible. In some portions of the 
galaxy where both high and low crcomponents are seen the 
extinction in the starburst component is slightly stronger than 
in the shocked component, though the average difference in 
E(B — V) is comparable to the intrinsic spread in values. 

After applying the velocity dispersion cutoffs, the total 
Ha flux of the high-er emission line gas accounts for as much 
as 1/3 of the total Ha emission, which would decrease an es- 
timated star formation rate by ~ 30%. Clearly shocked emis- 
sion must be considered when SFR measurements are made 
from Ha in galaxies with starburst driven winds and merger 
widespread shock excitation 

7.4. Aperture & Redshift 

As shocks contribute fractionally to the emission line fluxes 
over the majority of the galaxy and dominate the total flux in 
some off-nuclear regions, care must be taken when deriving 
quantities from the emission line spectra. In order to better 
understand the effect of aperture size and it's changing physi- 
cal scale with redshift, we artificially redshifted our data cubes 
while extracting a psuedo-fiber spectrum from the peak in me- 
dian rest-frame red light flux. 

We choose SDSS as our point of comparison, using a 
3" fiber size for our aperture extractions and rebinning our 
spectra in wavelength to match the approximate R~ 1800 
spectroscopic resolution of SDSS. We maintained a spaxel 
size of 1 "while spatially rebinning and dimming the data 
cubes to correspond to z of 0.2, 0.4, 0.8, 1.5 and 2.5. After 
extracting a composite pseudo-fiber spectrum at each redshift, 
we apply our fitting routine using both one and two compo- 
nent Gaussian fits to emission lines. The measured emission 
line ratios for each redshift and each component are shown 
in figure QT| Asymmetries in the emission lines are in some 
cases fit somewhat better by 2 components, but the overall 
separation between the two components is ~ 100 km s -1 , 
within the simulated resolution. 

As we redshift our data, the emission line ratios do not 
change by more than 0.4 dex at the most. Most of this shift 
occurs after the first redshift owing to the binning of a large 
part of the galaxies emission. Both galaxies move from a 
more composite to a more starburst-like classification in the 
[NII]diagram. The two component fit is more difficult to in- 
terpret, but generally the broader component in both systems 
moves to the right in the diagnostic diagrams (more shock 
like) while the narrower component moves to the left (more 
HII region like). Because of the low separation in velocity 
space the interpretation of the secondary component is proba- 
bly unreliable. 

With such an SDSS fiber aperture at the highest simulated 
redshifts (z=1.5, 2.5) the spectra become a sum of nearly all 
of the emission line gas in the system and the total emission 
line fluxes of the weaker lines decrease to a few x 10 _17 ergs 
s _1 cm~ 2 or less, making them marginally detectable. Pro- 
portionally more of the high-cr gas is of lower surface bright- 
ness, so that while more of the outer regions dominated by 



high-cr shocked gas are coadded, their total surface brightness 
drops, countering this. The total velocity dispersions of the 
narrow portion of the emission lines also becomes broader as 
more HII region ensembles are coadded, increasing to higher- 
a of ~ 100 consistent with those measured in high -z galaxies. 
This e ffect was noted using IFU observations by Epin at et al.l 
d20T0 ) and causes difficulty in discerning between large dis- 
persion clumps and effects and the coadding of several low-cr 
regions at high redshift. This is especially noticeable in IC 
1623 which has a complex velocity structure across the ma- 
jority of the brightest portions of the galaxy. 

Because of these effects, separating the shocked gas com- 
ponent from the star forming component within global spectra 
is difficult at the highest z. Improvement can be achieved with 
higher spatial and spectral sampling coupled with deep obser- 
vations, though the purely shocked high-a regions associated 
with the dust lanes in NGC 3256 and IC 1623 are of lower 
surface brightness making them difficult to detect even if the 
global spectra show composite behavior. 

8. CONCLUSION 

Using IFU emission line data of two nearby LIRGs, we 
have found a bimodal distribution of velocity dispersions and 
emission line ratios indicative of composite activity. This 
composite emission can be explained by a combination of 
shocks and star formation. In order to separate the shocked 
component from the purely star forming component we estab- 
lish cutoffs in velocity dispersion to distinguish the shocked 
higher-er component from the HII region-like lower-cr gas. 
Once this distinction has been made, the HII region like gas 
can be used to calculate derived quantities based on HII region 
emission line data such as star formation rate and metallicity. 
The properties and overall contribution of the shocked gas can 
also be inferred from the emission line data as well. We have 
provided such an example using low-velocity shock models 
which accurately reproduce the observed emission line ratios 
at the extreme end of the emission line ratio diagnostics. The 
shape of the mixing sequences between these models is di- 
rectly influenced by the gas abundance and our HH-region 
models imply a distribution of Q(H) consistent with with th e 
dynamical models for HII regions bv lDopita etail (I2006albl) . 

We urge caution when interpreting emission line data of 
both nearby and high-z galaxies, especially when spatial in- 
formation is limited. Traditionally the classification of a 
galaxy relies on a nuclear spectrum, or a globally extracted 
aperture in the case of more distant systems. Using this 
method our systems would be classified as composite locally 
and starburst at higher-z, though this classification is affected 
by the size and position of the aperture extracted and the mor- 
phology of the system. Slit spectra provide a slightly more 
complete picture though the alignment and size of the slit and 
the size and distance of the system may complicate interpreta- 
tion. This is especially true at high redshift (z~ 1.5 — 3) where 
a slight misalignment of the slit coupled with the drastic re- 
duction in surface brightness of most features makes such 
measurements difficult and may miss the complete picture. 

Our data clearly illustrates the importance of shocks in 
merging galaxies-not only do they affect quantities derived 
from spectroscopy, they also act as a useful mechanism for 
dissipating the kinetic energy and angular momentum of the 
infalling gas in merging systems. This aids in driving the gas 
to the center of the galaxy in order fuel the ongoing central 
starburst and keeping it in an increasingly deep potential well 
and presumably feed the eventual AGN. Presumably a por- 
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FIG. 1 1. — Single component fit to artificially redshifted, SDSS fiber aperture-extracted data cubes. Both IC 1623 and NGC 3256 are shown and points are 
color coded according to redshift. The points do not shift more than ~ 0.4 dex. NGC 3256 moves from a composite classification to an HII classification in the 
[Nil] diagnostic. 



tion of the shock excitation is generated by galactic winds 
as the merger progresses, though even a few percent of the 
total shock luminosity we calculate is sufficient to dissipate 
the binding energy on a time scale shorter than the dynamical 
timescale (10 to 100 Myr vs. 800 Myr). As the gas is con- 
tinually driven inwards and shocked by tidal forces, galactic 
outflows coupled with star formation will also compete to use 
up the gas on the scale of a few hundred million years or so. 

Our work also illustrates the utility of Integral Field Spec- 
troscopy in accurately decomposing the complex aspects of 
a merger. Not only is the full spatial picture of many prop- 
erties available with an IFU observation, it is easier to sepa- 
rate the widespread shocked component and to trace its spatial 
and kinematic structure. Once this information has been iso- 
lated, the derived quantities of the HII region emission can 
be more accurately calculated and the detailed properties of 
the shocked gas itself can be investigated. This is especially 
important when calculating metallicities from strong-line cal- 
ibrations, a quantity of particular interest in merging galaxies 
locally and in the early universe. A forthcoming paper (Kew- 
ley et al. in prep) will provide a theoretical framework for 
the effect of shock emission on such measurements and dis- 
cuss the mixing sequence in galaxies such as those presented 



in this paper. We also plan to apply this careful separation of 
shocked gas from HII region emission where necessary to the 
measurement of metallicities in a larger sample of IFU data 
of merging galaxies in a separate paper (Rich et al. in prep). 
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